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Abstract Cosmological numerical simulations of galaxy formation have led to the 
, cuspy density profile of a pure cold dark matter halo toward the center, which is 

in sharp contradiction with the observations of the rotation curves of cold dark 
matter-dominated dwarf and low surface brightness disk galaxies, with the latter 
■ tending to favor mass profiles with a flat central core. Many efforts have been 

devoted to resolve this cusp-core problem in recent years, among them, baryon- 
^ . cold dark matter interactions are considered to be the main physical mechanisms 

' erasing the cold dark matter (CDM) cusp into a flat core in the centers of all CDM 

, halos. Clearly, baryon-cold dark matter interactions are not customized only for 

CDM-dominated disk galaxies, but for all types, including giant ellipticals. We flrst 
flt the most recent high resolution observations of rotation curves with the Burkert 
CO ' proflle, then use the constrained core size-halo mass relation to calculate the lensing 

^ [ frequency, and compare the predicted results with strong lensing observations. 

Unfortunately, it turns out that the core size constrained from rotation curves 
' of disk galaxies cannot be extrapolated to giant ellipticals. We conclude that, in 

, the standard cosmological paradigm, baryon-cold dark matter interactions are not 

• ' universal mechanisms for galaxy formation, and therefore, they cannot be true 

I solutions to the cusp-core problem. 

I Key words: cosmology: dark matter — galaxies: formation — galaxies: halos 

galaxies: structure — gravitational lensing 



1 INTRODUCTION 

In the standard A cold dark matter cosmology paradigm (called LCDM), the observed large 
scale structure (LSS) of the universe, which is composed of virialized galaxies and clusters 
of ga laxies, is hierar chically formed from primordial density fluctuations in the early universe 
(e.g.. iLongaiii [20081 ). During the matter dominated epoch, CDM over densities provide deep 
enough gravitational potentials so that small dark halos can form, merge and evolve more 
rapidly (compared with only baryons) to make the present LSS. It is believed that, during 
the whole structure formation process, baryons follow CDM potentials and then reside at the 
centers of their host dark halos. As a result, each galaxy or cluster of galaxies is surrounded by 
a dark matter halo. In the framework of LCDM N-body numerical simulations, some analytical 
and semi-analytical models have successfully predicted most of the observational properties 
of LSS. However, challenges arise from small, galactic scales. Numerical simulations predict 
too many substructures (dwarf galaxies), e.g., several hundreds of satellites are predicted in a 
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galaxy like our Milky Way, while only 23 of them are observed; to o low angular momenta of 
spira l galaxies; and too high densities (cusps) in the galactic ce i iters ( Navarro et al.l . Il997l : iJind . 
I2OOOI: iNavarro et all . l2004l : Ijing fc Sutol . l2002t iGraham et all . l2006l ). whereas high resolution 
observations of the rotation curves for CDM-dominated dwarf and low surface brightness (LSB) 
galaxies imply that galactic dark matter halos have a density profile with a flat central core 
(Ide Blok et al.l. l2001bl: Ivan den Bosch fc Swate^. I2OOII: ISw atcrs e t all. l2003t IWeldrake et all 
I2OO3I: iDonato et all . 12004 iGentile et al.l . l2005t ISimon et al.. . .20051 iGentile et all l2007l ). The 
latter is known as the cusp-core problem, which we will investigate in detail in this paper. 

Many efforts have been devoted to resolve the challenges mentioned above. The less tra- 
ditional models include warm dark matter and self-interacting dark matter. In the LCDM 
model, some researchers simply deny the validity of the observations, others introduce 
baryon-cold dark matter interactions (BCDMIs) such as dynamical friction of substructures 



(lEl-Zant et al.l. I2OOII: iToiiini et all. l2006t iRomano-Diaz et al 
t ionj Weinberg fc Katz j_ 20021: Hollev-Bockelmann et al 



(Mashchenko et al. 



200a 



iPeirani et al. 



,120081). 



2008,), stellar bar-CDM intcrac- 
2005f l. and baryon energy feedback 
In this paper, we focus on the validity of 
BCDMI solutions to the cusp-core problem. 

First of all, we point out that any reasonable solutions to the cusp-corc problem should 
be verified by the observations of galaxies not only on small scales like dwarf and LSBs, but 
rather on all galactic scales, especially on scales like giant ellipticals. Why should we include 
the giant ellipticals when testing structure formation models (e.g., the solutions to the cusp- 
core problem)? According to the hierarchical structure formation theory in the frame work of 
LCDM, small CDM halos form first, then merge to make increasingly large halos. While N-body 
numerical simulations of pure CDM with higher and higher force and mass resolution still favor 
cuspy halo profiles, the BCDMIs introduced in the baryon plus CDM regime can transfer the 
early formed cusp into a central flat core in later stages for each dark h alo. It is well known th at 
the "universal" density profile known as NFW (Navarro- Frenk- White; iNavarro et al.l . Il997t) is 
valid for halos ranging from dwarf galaxies to those as large as clusters of galaxies. The only 
difference between pure CDM and baryon plus CDM regimes is that the latter considers the 
effects of baryons. Their formation processes should be similar (bottom-up or hierarchical) in 
the spirit of LCDM. Therefore, the properties of the galaxies (e.g., the slope in the central 
region) predicted by BCDMIs should be valid and meaningful at all mass scales. Needless to 
say, for any theoretical models describing structure formation, it is meaningless or impossible 
to customize some physical processes only to explain the observations of dwarfs and LSBs, and 
exclude ellipticals as exceptions. In other words, no matter what galaxy types the initial con- 
ditions presumed, the subsequential and final (present) CDM halos predicted by simulations 
should include all galaxy types and have the same observational properties. This is why we 
should include giant ellipticals when testing solutions to the cusp-core problem. Usually, when 
we talk about the cusp-core problem, we mean the contradiction between the "cusp" predicted 
by simulations of pure CDM and the "core" implied by rotation curves of dwarfs and LSBs, as 
is widely known. However, we forget the fact that giant elliptical galaxies are also important 
members of the galactic family, but they prefer a cuspy (or small core-size) density profile in 
their central region according to observations (e.g., x-rays and the statistics of strong gravi- 
tational lensing). In deed, LSB galaxies and giant ellipticals are quite different galaxies; the 
former are CDM-dominated while the latter have baryon-dominated centers. However, their 
distinguishing differences are obtained from all kinds of observational data rather than theoret- 
ically predicted from structure formation theories. It is theories that should successfully explain 
such observational differences, not the other way around. The density profiles of giant ellipticals 
with baryon-dominated centers (a cusp or small core) should be predicted by the same model, 
if reasonable, which predicts the density profile of LSB galaxies with CDM-dominated centers 
(a large flat core). 

It is interesting to note that some solutions to the cusp-core problem based on analytical 
models are problematic in a similar way to simulations. One such model assumes that cosmo- 
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logical halos form from the collapse and virialization of "top-hat" density perturbations and 
are spherical, isotropic and isothermal. This predicts a unique, nonsingular, truncated isother- 
mal sphere (NTIS) and pro vides a simple phys i cal clue about the ex istence of soft cores in 
halos of cosmological origin (IShapiro et al.l . ll999l : lliie"v fc Shapiro! . l200lD . It is claimed that this 
NTIS model is in good agreement with observations of the internal structure of dark matter- 
dominated halos on scales ranging from dwarf galaxies to X-ray cluster s. Unfortun ately, the 
NTIS model is ruled out by the statistics of strong gravitational lensing |Che3,[200i). 

In summary, any self-consistent, reasonable solutions to the cusp-core problem, either nu- 
merical or analytical, should be able to successfully explain the cuspy density profiles of giant 
ellipticals simultaneously, otherwise, no matter how reasonable it looks, such a solution cannot 
be true. 

While giant ellipticals should be included, now the question is: to what extent is the core 
size of the giant ellipticals allowed by observations? In particular, is it acceptable to extrapolate 
the core size-halo ma ss relation constrained by t he rotation curves of dwarfs and LSB galaxies to 
giant ellipticals? Like lChen fc McGaugh I ( 2008[ ). we investigate these problems by "reasonably" 
extrapolating the core-mass relation constrained from LSB galaxies to strong lensing galaxies 
(usually ellipticals), and compar ing the predicted lensing frequency with observations. Since 
it was explicitly pointed out in iMashchenko et al~l ( 2006[ ) that the CDM density profile of 
the evolved halo is in remarkable agreement with the Burkert profile, we fit the observational 
data of rotation curves with the Burkert density profile in section [21 the corresponding lensing 
probabilities are presented in section [31 and we give our discussion and conclusions in section 

SI 



2 DENSITY PROFILE CONSTRAINED FROM OBSERVATIONS OF 
ROTATION CURVES 



We use high- resolution a nd high-quality hybri d Hq/H I rotation curves of a sample of 26 
analyzci 

profile (jBurkerO . [l995t ) . S ome curves of th i s sam ple were taken from the large sample of 50 



We use nign- resolution a nd nign-quality nybri d -tiQ,/Jrl i rotation curves oi a sample or zb 
LSB ga laxies, an alyzed in' de Blok et al.l ( 2001a! ). to fit the parameters in the Burkert density 



LSB galaxies presented in lMcGaugh et all (!200ll). While these curves were well fitted by the 
cored isothermal sphere model (jBegeman et al. . 19911 : !de BloK fc Bosmal [2002t !de Narav et all 
!2008l). we will show bellow that they are also well fitted by the Bur kert density profile (see 
also ISalucci fc Burkertl . !2000l: iGentile et all . 12004! : ISalucci et all . !2007t) . As a typical model of 
BCDMIs, turbulence d riven by stellar feedback du ring galaxy formation is a possible solution to 
the cusp-core problem ( Mashchenko et al~l . !2006[ ). Numerical N-body simulations show us that 
random bulk motions of gas in small primordial galaxies result in a flattening of the central CDM 
cusp. Phase space arguments imply th at the core should persist through subequent mergers 
( Dehnenl . 120051 : [Xazantzidis et al.L!2006f ). Consequently, in the present universe both small and 
large galaxies would have fiat CDM core density profil es, which can be well fitted with the 
Burkert model (|Burkertl . [19951: IMashchenko et al. I . l2006l) . 



p{r) 



(r + ro) (r2 + r^) 



(1) 



where po is the central density near r = and rp is the cor e radius, the two free parameters to 
be fitted. The mass within radius r is ( Salucci et al.l . [20071 ). 



M{r) = 2TTparl 



In ( 1 + — 



In 



1+ ( - 

^■0 



— arctan ( — 




(2) 



The corresponding orbital velo city at radius r is simp ly V'^{r) = GM{r)/r, where G is the 



npi; 

gravitational constant. The sample (|de Blok et al.l . l2001a[) provides for each galaxy the observed 
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velocity y/(Km/s) and uncertainty yErr/(Km/s) at radius r/Kpc. We fit the rotation curves 
with V^{r, po,rQ) and derive the values of the parameter pair {po, rg) for each galaxy of the 
sample. Figure [T] shows the fitted V'^{r) (solid lines) and the observed data of rotation curves 
(crosses) for each galaxy. Note that, for a few galaxies, especially F579, U11648 and U11748, 
the Burkert-fitted lines (solid) deviate greatly from the data points at large radii. However, this 




Fig. 1 Comparison of the fitted V'^{r) derived from the Burkert profile (solid lines) 
and the observed rotation curves (crosses). 
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Fig. 1 Continued. 



will not affect our further analysis, because our interests are focused on the central regions, 
where the data have been well-fitted. Statistically, there is a correlation between po and vq for 
the sample; the fitted result is 

X -1/2 
^0 \ , ^3 



Po= 0.029(^) M0/pc^ (3) 



In later lens probability caculations, we need to know the total virializcd halo mass M ~ My 
For the Burkert profile, the mass diverges logarithmically at large radii, and should be cut off 
at r200: the radius of a sphere within which the average m ass density is 200 tirn es the critical 
density of the universe, typically taken as the virial radius ( Navarro et al.l . [l997t ). So we define 



,,r ,,r/ ^ 8007rpcrit(z) 3 

M200 = M{r2m) = ^ ?-2oo. (4) 

where M(r2oo) is the virialized halo mass obtained by substituting r = r2oo into equation 
and Pcrit(-z) is the critical density of the universe at redshift z. By numerically solving Equations 
^ and (H]), we get an approximate relation between rg and M, 

M 




Fig. 2 Left panel: the correlation between central density po and core radius vq; crosses 
represent the data derived from the rotation curves of the sample with the Burkert 
profile, the solid line is the fit with the data. Right panel: the relation between core 
radius rp and the virialized halo mass M . 



which is similar to the result of ISalucci et al. 1 (l2007f ). although we have adopted different pro- 
cedures and used different samples. The correlations between the pairs (po, tq) and (tq, M) are 
presented in Figure [21 



3 STRONG LENSING PROBABILITIES 



It has been explicitly pointed out ( Mashchenko et al~l . l2006l ) that, on one hand, the resultant 
cored dark matter halos produced by stellar feedback are in agreement with the Burkert profile, 
and on the other hand, as the dwarf galaxies merge together to make larger ones, the flat-cored 
shape of the dark matter density profile is preserved in subsequent massive halos. Therefore, 
in this section, we assume that the Burkert profile constrained from the rotation curves of the 
LSB galaxies in the previous section is also appropriate for strong lensing galaxies (usually giant 
ellipticals), and we use it to calculate the corresponding lensing frequency. 

Gravitational lensing p rovides a power f ul and independent tool t o det ec t th e mat- 



ter d istrib u tion in LS S (ISchneider et al" . 1992^: ^Wambs ganss et ahl. 
20001 12004 iBartelman^. Il996l: iBartelmann fc Schn eider. 200 1|: iKeeto"!]! . 



1995 



2001 



1996 



2003 



Keeton fc Zabludofi . 2004 Zhang et al. . 2003 . |2009[ ). In particular, by virtue of eliminating 



the degeneracy of some lens and cosmological mo dels, strong lensing statis t ics exhibits an ex- 



ceptional ability to invest i gate both dark matter ( Turner et al 

2002; 'Li fc Ostrike^. |2002|. [2003 IChad. [2005 IChe n. '2003a"b' 200. , . .. . . 

200& .Qguri et all. I2002|l0gurll2003t lOguri et al.,. .2003; .Zhang . ,2004 IChad. l2007l) and dark 



1984 IWul. 119931: IChae et al 
' IChen fc Ziiao. l2006l: |C~ 



energv (|Fukugita et all . Il990l : iFukugita. fc Turneil . ll99lHTurneI Il990t iKrauss fc Whit3 . Il992l : 
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19991: IWaga fc Micelil. Il999l: ISarbu et"^ 



1995 



2001 



1996HFalco et all. Il998l: ICoorav fc: Huterer 
Chenl . l2004allbl : iDev. Jain k Alcanid . 12004 



Yang fc Chenl. l2009HZhang et all . 12009^ With the increasing number of the newly discovered 



lenses (e.g., Wen et al.l . 20C)9() . we now have at least t wo weh-defined statistical lens samples . 
In addition t o the radio lens sample (C LASS/ JVAS ; iPatnaik et all . Il992l : iKing et all . Il999l : 
iMvers et al. ], l2003l : iBrowne et all . l2003l). we have a Sloan Digital Sky Survey Quasar Lens 
Search Data Release 3 (SQLS DR3; Inada et al.l . 120081 ) optical lens sample, which first included 
a cluster-scale lens, and thus provides stronger constraints on halo models, espe cially for the 
universal properties of halos ranging from dwarf galaxies to clusters of galaxies ( Oguri et al.l . 
[2008; Oguri & Blandford, 2009). 

The lens equation is 77 = -DsC/^'l — ^lsQ^, where rj and ^ are the physical positions of a 
source in the source plane and an image in the image plane, respectively, a is the deflection angle, 
and Z^L, ^^id Z^ls s-re the angular diameter distances from observer to lens, observer to 
source, and lens to source, respectively. By defining dimensionless positions y = Di^rj / DgrQ and 
X = ^/rg, and dimensionless angle a = Di^Di^ga/ Dgro, the lens equation is then y = x — a. For 
circularly symmetric le ns, a = AGM{^) / c^^^ where M(^) is the mass within a circle of radius £ 
(jSchneider et al.l . ll992h . For the Burkert profile, M{x) = ^■nporlm{x) (jPark fc FergusonLl2003l ). 
where 



m{x) 



In f + ^(VP TT - 1) + ^^ ^arccoth(V?TT) 
— i ^/ x"^ — 1 arctan — 1, 

-ln2-f + ^[7r + ln(3 + 2V2)], 
In f f (VPTT - 1) + ^^iSIarccoth(^/Pn) 



if a; > 1 



if; 



1, 



- ^ Vl — x^arctanhVl — 2;^, 
The lens equation for the Burkert profile then is 



if X < 1. 



y = x 



Ikc m{x) 



27r^Gporo \ f DiJ^\ ^ S(0) 



J \ Ds J 



■^crit 



(6) 



where Kc is the central convergence, an important quantity that determines wheither strong 
lensing can occur or not, S(0) = 7rpo?'o/2 is the surface mass density at the lens center and 
^crit ~ Ds/ 4:TtGDi^Di^s is the critical surface density. 

Generally, for any spherically symmetric density profiles of lensing hal os, multiple image s 
can be produced only if the central convergence Kc is greater than unity ( Schneider et al.l Il992[ ) . 
When Kc < 1, only one image is produced. Note that even if Kc > 1 is satisfied, mult iple images 
can occur only when the source is located within j/cr = y{xcr) l|Li fc Ostrikerl . [2002^ . where x^r 



wi th dy/ dx = for a; < (this is similar to lensing by cored 
Che"3 . 120051) ■ For a singular density profile such as the singular 



is determined from equation 
isothermal sphere halos, see 
isothermal sphere (SIS) and NFW profiles, the central value is divergent, so k > 1 is always 
satisfied, and multiple images can be produced for any given mass. For density profiles with a 
finite soft core such as core isothermal sphere, NTIS and Burkert profiles, however, the condition 
K > 1 requires that only halos with mass greater than a certain value (determined by Kc = 1) 
can produce multiple images. 

The calculations for lensing probabilities with quasars at redshift lensed by Burkert h alos 
of galaxies and with image separations larger than A9 are straightforward (see IChenl . l2005l . for 
details) , 



P(> A0) = OMyc 



(1 + 4 



,dD£izl 
dz 



dz 



f{M, z)dM 



y[^i{y)Y■'dy, (7) 
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CLASS/JVAS 
SIS in LCDM 
BURKERT in LCDM 
NFW in LCDM 




0.5 1 1.5 2 

A0(arc seconds) 



2.5 



Fig. 3 Lensing probability with image separation larger than A0. Our predicted 
lensing probability based on the Burkert model (dotted line) is about four orders 
of magnitude lower than the observations of CLASS/JVAS (thick hostogram). As 
expected, the SIS model (dashed line) is in veru good agreement with observations. 
The NFW model (dash-dotted) also fails to explain observations but is much better 
than the Burkert model. 



where D£(z) is the proper distance from the observer to the l ens located at redshi ft z; f(M, z) is 
the mass function, for which we use the expression given bv I Jenkins et al.l ( 200l[ ): M^i^(A9, z) 
is the minimum mass of halos above which lenses can produce images with separations greater 
than A0; and ^(y) is the total magnification of the two outer images for a source at y, which 
can be calculated directly from the lens equation 

In our actual numerical calculat ions, we cho ose almost the same procedures, mass function 
and even background universe as in Chen! (l2005ll, except for the distribution of the redshifts of 
source quasars, for which, as in lChen fc McGaueh I ^M), we use only the mean value Zs = 1.27. 
The result is shown in Figure [31 The observational result for the well defined combined sample 
of the cosmic Lens All-Sky Survey (CLASS) an d Jodrell Ba nk/Very Large Array Astrometric 
Survey (JVAS) is shown as thick histogram fsee lCheiil . l2005l . and references therein). For com- 
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parison, the lensing probability of the SIS model and NFW model are also shown with the same 
parameters and approximations as Burkert. 



4 DISCUSSION AND CONCLUSIONS 



We have calculated the lensing probability based on the Burkert profile model, for which the 
core size-halo mass relation is obtained by fitting the high-resolution rotation curves of a sample 
containing 26 galaxies. For comparison, we have recalculated the lensing probability for the 
SIS model and NFW model with the same parameters and approximations as Burkert. As 
expected, the lensing probability for the Burkert model is 4 orders of magnitude lower than the 
observations of CLAS/JVAS. While our stong lensing calculations rule out the Burkert model as 
a possible gia nt elliptical density profile , we conclude, at the s ame time, that the sellar feedback 
mechanisms ([Mashchenko et aL 1. 120061 : If eirani et al. 1. 12008D are also ruled out as reasonable 
solutions to the cusp-core problem. Additionally this conclusion is independe nt of the models 
used to fit rotation curves: Burkert in this paper, and cored isothermal sphere ( de Narav et al.l 
I2008D in lChen fc McGaueh I ^QOi)- 

Other solutions to the cusp-core p r oblem , incl uding the analytical NTIS model 



( Shapiro et al 



(El-Zant et al 



19991 : Eiev fc Shapirol. 20011: Chenl. 2005), dynamica l friction of substructures 



20011: iTonini et all. 20061: iRomano-Diaz et al. I. 20081) and stellar bar-CDM in 



teraction ([Weinberg fc Katz I. l2002r iHoUev-Bockelmann et al.H . |2005| ). face exactly the same 
embarrassing problem. Clearly, the solutions proposed so far can produce a large core for each 
dwarf and LSB galaxy, and thus can successfully explain the observations of rotation curves, 
but they cannot explain the steep and cuspy centers of massive galaxies, which are favored by 
stong lensing and X-ray observations. As a matter of fact, in the framework of LCDM, the so 
called "cusp-core problem" encountered by hierarchical structure formation theories is not just 
the contradiction between the cusp predicted by pure CDM simulations and the core implied 
by the observations of rotation curves, but rather, the cusp predicted by pure CDM simulations 
on all mass scales, the core implied by rotation curves on small, dwarf galaxy scales and the 
cusp favored by strong lensing and X-ray observations on large, giant elliptical galaxy scales. 
From an observational point of view, regardless of theories, we can call this the "cusp-core phe- 
nomena" rather than "problem" , when referring to the cusp in large scale galaxies and core in 
small scale galaxies. As we have pointed out, a reasonable, self-consistent theory should explain 
the observations on all mass scales, i.e., all "cusp-core phenomena", rather than only the cores 
on small scales. 
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